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This work presents an upper bound on the neutrino mass using the emission of v e from the 
neutronization burst of a core collapsing supernova at 10 kpc of distance and a progenitor star of 
15 Mg. The calculations were done considering a 34 kton Liquid Argon Time Projection Chamber 
similar to the Far Detector proposal of the Long Baseline Neutrino Experiment (LBNE). We have 
performed a Monte Carlo simulation for the number of events integrated in 5 ms bins. Our results are 
< 2.71 eV and 0.18 eV < m „ < 1.70 eV, at 95% C.L, assuming normal hierarchy and inverted 
hierarchy, respectively. We have analysed different configurations for the detector performance 
resulting in neutrino mass bound of 0( 1) eV. 


I. INTRODUCTION 

From the current experimental evidence of neutrino oscillations mm, we know that neutrinos are massive particles, 
and each neutrino flavor eigenstate \vp) (/? = e, \x, r) is a coherent superposition of mass eingenstates | v.j) (i = 1, 2,3), 
connected by the PMNS mixing matrix Upi, through the relation \vp) = JT Upt\vi). However, oscillation experiments 
are sensitive only to the mass squared differences (A = m 2 — m 2 ; i,j = 1,2,3) and not on the absolute 
value of each mass eigenstate which is one of the most challenging measurement in particle physics. Furthermore, 
one must be careful about the interpretation of experiments on neutrino masses since different techniques result 
in different observables. For instance, neutrino-less double beta decay experiments search for Majorana neutrinos 
with mass m 2 e = \Y2iU 2 i m i\ 2 an d single-/? decay experiments observe features driven by massive neutrinos with 
m 2 = \U 2 i\ 2 rn 2 not dependent on Majorana phases. Since m e > m ee measurements which are not exclusively 

sensitive to Majorana neutrinos can be safely taken as conservative ones [3]. 

If we can measure the mass of any neutrino flavor up the masses of v,j can be obtained from the boundaries of the 
mixing matrix and other A m 2 j obtained from different oscillation experiments. 

The strongest experimental limit in the neutrino mass, as quoted in the Particle Data Group (PDG) [4], was 
obtained by the Troitsk experiment on tritium single-/? decay: m e < 2.05 eV at 95% confidence level (C.L) [5]. This 
kind of experiment measures the modification of the /? spectrum near its endpoint caused by the neutrino mass. The 
KATRIN experiment 6] will improve the technique and plan to set a stronger limit by enhancing the sensitivity at 
least in one order of magnitude. 

Neutrinos emitted from stellar collapse, for example type II supernovae, can give a valuable contribution to the 
determination of the mass of neutrinos, as idealized by Zatsepin [7]. Using neutrino data [8HTU] from the SN1987A, 
which exploded in the Large Magellanic cloud (~50 kpc), and the likelihood of detection analysis proposed in [III , 
an upper bound of 5.8 eV at 95% C.L m was obtained. With a different likelihood approach and astrophysical 
parametrization of the neutrino emission, Lamb & Loredo m reached a 5.7 eV at 95% C.L for the neutrino mass 
limit. 

There are good perspectives to explore the sub-eV region in the future experimental scenario for detection of galactic 
supernovae muMi]. Other mass limits obtained from the supernovae can be seen in the table presented in Ref. |18| . 

The perspective to detect v e from a supernova, but in a liquid scintillator detector, was recently considered in using. 
Water Cherenkov detectors with Gadolinium could also improve the number of u e events according to [2T| . In the 
present work, we aim to discuss bounds on neutrino masses from supernovae motivated by the Zatsepin’s idea. We 
are going to consider electronic neutrinos (u e ) emitted from the neutronization burst of a galactic supernova at a 
distance of 10 kpc and 15 M e and detected by a Liquid Argon Time Projection Chamber (LArTPC). 

The paper is organized as follow: Sec. |n] presents the main assumptions about electronic neutrino production and 
its propagation; Sec. [Hd] presents the basic features about the neutrino detection in the LArTPC; Sec. [TV] presents 
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the main assumptions about our event generated Monte Carlo simulation and Sec. [V] presents the boundaries on the 
neutrino mass. Sec. VI and Sec. VII present the discussion about our results and our main conclusions, respectively. 


II. NEUTRINOS FROM SUPERNOVAE 

As the shock evolves against the infalling dense matter of the outer core of the star, some of its energy is lost 
by the photodissociation of nuclei into protons and neutrons. At this moment, there is an abundant production of 
electronic neutrinos (u e ) via p + e~ —> n + u e . These electronic neutrinos accumulate in a very dense and opaque 
environment behind the shock wave. When the shock reaches the zone with densities around 10 n g cnU 3 in a few 
miliseconds after the bounce the i/ e ’s are released. This process is called neutronization burst. It has a peak luminosity 
of ~ 3.5 x 10 53 erg s _1 and lasts ~ 25 miliseconds. For a recent review of the current knowledge and status of the 
explosion mechanism, see [ 22 ] ■ 

At a distance d of the progenitor star, for each flavor up (/3 = e, e, x), the unoscillated sprectral number fluxes 

( F °p(E)) are: 


170 ( T?\ __ A,, / "(I (E) 


( 1 ) 


where L Vfj is the luminosity for the respective flavor up, (E Vfj ) is the up mean energy and f Vff (E) is the quasi-thermal 
spectrum written as follows B5| : 

Up{E) = X/3 (^y Y e -(^ + i)E/(E^)_ (2) 

The ap is a parameter defined by {E% ) /(F ^) 2 = (2 + ap)/{l + ap) and \p is the normalization constant factor: 
fdEf Vp (E) = 1. Several of these parameters, such as L Vfj in Eq. Q, and (F 2 ^) and {E Vp ) 2 in Eq. Q, change their 
respective values with time after bounce and depend on the simulation of the star explosion [24j . 

After taking into account the effect of neutrino oscillations, in the three neutrino families framework in a supernova 
environment, the fluxes (F Vfi ) at Earth for each flavor up are [251 : 


F Ve =pFl + (l-p)F^, 
Fp e =pFl + (l-p)F? x 


(3) 

(4) 


4F„ x = (1 - p)F° + (1 - p)Fl + (2 +p + p)F ° x , (5) 

where are the primary unoscillated neutrino fluxes at the production region, x = p or r and p(p) is the survival 
probability of an electron (anti)neutrino after propagation through the SN mantle and the interestellar medium. For 
the actual oscillated parameters, p = 0 for normal hierarchy (NH) (toi < m 2 <C m 3 ) and p = sin 2 $12 for inverted 
hierarchy (IH) (m 3 <C m 2 < mi). 

For massive neutrinos, the arrival time (t®) of the neutrino in the detector at Earth will be delayed [7] : 

i© = ti + d/c + At, (6) 


where ti is the emission time at the source and c is the speed of light. The delay introduced by the neutrino mass m v 
in Eq. ([d]) is represented by At and can be written as: 




(7) 


III. DETECTION ASSUMPTIONS 

The neutrino flux of the neutronization burst is mainly populated by electronic neutrinos. Therefore, only experi¬ 
mental technniques with high sensitivity to u e can have a possibility to investigate this particular phase of supernova 
explosion. 
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The Liquid Argon Time Projection Chamber (LArTPC), such as the ICARUS-like detectors [55], have all the 
features required to identify the neutronization burst. The “Long Baseline Neutrino Experiment” (LBNE) [27: 
foresees a 34 kton (fidutial volume) LArTPC as its far detector. The experiment is dedicated for determination of CP 
violation phase and neutrino mass hierarchy. Though, its dimensions and underground site are suitable for exploring 
rare events, such as neutrinos from supernovae and proton decay. This experimental design is taken into account for 
our calculations. 

The most favourable electronic neutrinos ( v e ) detection channel in LArTPCs is the charged current (CC) interaction: 
L' e + i0 Ar —> e~ ± 40 K*. There are other detection channels, but all of them with smaller cross section. These detection 
channels can be tagged by spectral analysis of the photon emission from de-excitation of K , Cl or Ar , which exhibit 
specific spectral lines. The v — Ar cross sections can be found in Fig. 3 of Ref. [28. . 


IV. CALCULATION METHOD 

The number of expected events in the detector can be evaluated in the following way: 

J dEa(E) j dtiF Vfj (E,ti)5(tQ -U- At), (8) 

where cr{E) is the v e — Ar CC cross section, nt is the number of targets in the 34 kton Argon detector and F Vfi (E , ti) 
is the oscillated neutrino flux, which is evaluated by Eqs. (Jl]|5|. 

We performed a Monte Carlo simulation (MC) with very optimistic assumptions. Nevertheless, all of our assump¬ 
tions have solid grounds and can be taken as a “realistic case”. The number of expected events depends on the 
considerations we make about the neutrino production and their propagation till the detector. We considered a com¬ 
prehensive knowledge of the astrophysical process related to the SN explosion, which means that the astrophysical 
parameters related to neutrino production, ie the neutrino fluxes, are well known. However, to introduce uncertainties 
related to the explosion mechanism, we considered a 5% error in the overall normalization of the evaluated counting 
rates. Concerning the SN distance, we assumed a galactic SN of 10 kpc. The detector configuration, as pointed 
out before, is 34 kton LArTPC with 80% of efficiency for v e charged current (CC) detection. For now, background 
effects have been disconsidered since the detector is assumed to be placed underground. In this situation we expect a 
suppression of seven orders of magnitude [29]. The v e energy threshold for detection was set to 5 MeV. For the energy 
resolution we used a parametrization given by: 11 %/^E(MeV) + 2% [50] , For the time resolution, we adopted a 
±20 fis flat distribution. Also, we stress the fact that we performed our MC ignoring aditional complications given 
by the lack of knowledge if the first detected event is syncronized or not with the beginning of the burst. In the MC 
we included the two possible scenarions for mass hierarchy, as mentioned in Sec. [TTJ with the best-fit values of the 
oscillation parameters with no CP violation effects ED- 


V. RESULTS 


Our MC simulation results in thinner lines are shown in Fig. [TJ NH and IH hierarchy are represented by solid and 
dashed lines, respectively. The error bars in each bin are assumed to follow a Poissonian distribution and each bin 
has St = 0.005 s width. In Fig. [TJ we also show the expected number of events in the thicker lines considering a 0.5 eV 
neutrino mass, just for comparison with our MC. 

Non vanishing neutrino mass causes a delay in the experimental signal, changing the expected number of events in 
each time bin, as we can notice analysing Eq. (J?]) and Eq. ([8]). These theoretical modifications can be compared with 
the number of events in each bin generated by our MC simulation using the y 2 test: 


2 A (N t MC - Ni{m v )) 2 

X ~ 2-^ (T? 


(9) 


where i = 1, N bins, Nf AC is the number of events generated by our MC in the i-bin, Ni(rn u ) is the number of events 
evaluated by Eq. 0 and Eq. ([8]) for different neutrino masses {m„) and eq is the associated error. 

The “realistic case” provided Xmin ~ 2.12 for a neutrino mass of m„ = 0.00 eV and considering NH. From the 
Ay 2 = x 2 (w„) — Xmin > the upper bound on neutrino mass from the v e neutronization burst stage is, at 95% C.L, 
2.71 eV. On the other hand, for the IH case, we have obtained Xmin ~ 5.43 for a neutrino mass of m v « 1.00 eV. At 
95% C.L, 0.18 eV < rn v < 1.70 eV. We present our Ay 2 in terms of the neutrinos mass, m !/ (eV), in Fig. [2j for both 
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FIG. 1: (Colours online) Results of our Monte Carlo realization. For further details see text. Normal hierarchy (NH) in solid 
curves and inverted hierarchy (IH) in dashed curves. Thicker lines represent the expected number of events for a 0.5 eV neutrino 
mass. 
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m v (eV) 

FIG. 2: (Colours online) Neutrino mass bound considering a supernova at 10 kpc and a liquid argon detector with 34 kton, 
such as the one projected in the LBNE experiment. Dotted line is lcr limit; dashed line is the 2a limit; and dotted-dashed line 
is the 3 ct limit. The strong red dashed curve is for IH and the strong solid black curve is for NH. 

hierarchies. In these plots, dotted lines are bounds with 68.27% C.L, dashed lines are for 90% C.L; and dotted-dashed 
lines are the 95% C.L. 

We also performed other MC simulations and did the \ 2 analysis for the mass as a parameter for other situations 
that are described in the following: 

1 . “perfect case”: we consider that the normalization is completely known, 100% of detector efficiency and a perfect 
time and energy resolution. We mantain an energy cutoff in 5 MeV; 

2 . “perfect norm”: returning to the conditions estabilished in “realistic case”, but now for a completely known 
supernova normalization; 

3. “perfect eff”: considering a detector with 100% of efficiency and other conditions of the “realistic case” main¬ 
tained; 

4. “cut ener”: keeping the “realistic case”, but instead cutting events with energy less than 5 MeV, we perform a 
energy cut at ss 2.2 MeV; 

5. “time resol”: considering a temporal resolution of 1 ms instead of 20 /us of the “realistic case”; 

6 . “binl”: if we increase the bin size from 5 ms to 10 ms using the “realistic case” and 
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7. “bin2”: if we decrease the bin size from 5 ms to 2 ms using the other conditions of the “realistic case” 

We present in Table [T] our mass boundaries at 68.27%, 90% and 95% C.L for NH. The numbers in parentheses are 
the boundaries for IH. In this table we also show the minimum value of the mass ( rrimin ) with its respective y 2 , the 
y 2 ■ . 

A mm 



68.27% C.L 

90% C.L 

95% C.L 

'Wl'min 

Xmin 

realistic 

1.75(0.76-1.35) 

2.43(0.46-1.58) 

2.71(0.18-1.70) 

0.00 (1.00) 

2.12 (5.43) 

perfect norm 

0.91-3.25(0.65) 

3.79(0.95) 

4.07(1.1) 

2.27 (0.00) 

8.51 (3.22) 

perfect eff 

2.01(0.82) 

2.52(1.08) 

2.76(1.20) 

1.10 (0.00) 

2.39 (3.06) 

cut ener 

2.05(0.69) 

2.58(1.03) 

2.87(1.20) 

1.22 (0.00) 

7.53 (6.34) 

time resol 

1.34(0.71) 

1.94(1.03) 

2.22(1.17) 

0.00 (0.00) 

2.33 (8.53) 

binl 

0.71-2.23(0.63) 

2.66(0.93) 

2.88(1.09) 

1.54 (0.00) 

0.51 (3.22) 

bin2 

1.65(2.24-2.67) 

2.14(2.03-2.91) 

2.38(1.90-3.06) 

0.81 (2.50) 

4.84 (9.74) 

perfect 

3.26(0.96) 

3.74(1.19) 

3.96(1.29) 

2.37 (0.50) 

12.20 (4.41) 


TABLE I: Upper Bounds on m*,(eV) considering a supernova at 10 kpc and a liquid argon detector with 34 kton. We show the 
results for normal hierarchy (NH) and the numbers inside the parentheses are the mass evaluation for inverted hierarchy (IH). 
Also we show the the minimum value of the mass (m m i„) with its respective y 2 , the y ‘min- For details of the meaning of each 
test, see the text. 


VI. DISCUSSION 

From Table [i] we notice that, despite we change the configuration of the detector, statistical fluctuation of the 
Monte Carlo simulation plays a fundamental role in the mass bound determination, since we have few events in the 
neutronization burst and this process lasts only ss 25 ms. A very good example of this kind of behaviour is the one 
shown in the IH curve of Fig. |T| (red dashed). In the interval of 35-40 ms, there is a peak that was not theoreticaly 
predicted. So, this generates a worse y 2 and an interval in the neutrino mass bound (0.18-1.70 eV) with a Xmin 
1.00 eV. Another example of the impact of this statistical fluctuation is that, for the NH case, we have obtained a 
more stringent bound in the “realistic” case compared to the “perfect” case. 

Our bound, at 95% C.L, for the NH is very similar to the one obtained by the Troitsk experiment of tritium [3 
decay. So, if we really have NH and comparing with terrestrial experiments, which are independent of mass hierarchy, 
like Katrin (to*, < 0.2 eV), it will be very difficult that LBNE will put a more stringent limit on neutrino mass. 
From the Planck Cosmic Microwave Background measurements and galaxy clustering information from the Baryon 
Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey-Ill, and assuming a ACDM cosmological 
model, Ref. [32] obtained, at 95% C.L, ?7i„ < 0-35 eV. Ref. [33] obtained that Y2 m v < 0-18 eV, when one takes 
into account observations of the large-scale matter power spectrum from the WiggleZ Dark Energy Survey, Planck 
sattelite and baryon acoustic oscillation. If we consider that mass eigenstates are degenerated, m„ < 0.12 eV [32] and 
to*, < 0.06 eV (33]. 

For the NH and IH case, our bounds, respectively, at 95% C.L, are to*, < 2.71 eV and 0.18 eV < to„ < 1.70 eV. 
These bounds are almost similar with other prospects of neutrino mass bounds obtained for SN neutrinos. Just 
for comparison, authors in [34) considering only normal hierarchy, using a 20 kton liquid scintillator at Jiangmen 
Underground Neutrino Observatory (JUNO) and a galatic supernova v e signal obtained a limit on the absolute 
neutrino mass scale: m„ < (0.83 ± 0.24) eV at 95% C.L. This limit is also similar with the ones obtained in the 
other cases discussed in this work - see Table [TJ So, our conclusion here is that, for a SN at 10 kpc of distance and 
15 Mg, LBNE 34 kton argon detector can probe neutrino masses of order of 1 eV, at 95% C.L. We stress the fact 
that knowing neutrino mass hierarchy is fundamental to put limits on neutrino masses from SN and LBNE has as 
one of its objectives determine this property. 


VII. CONCLUSIONS 

Even though our supernova neutrino bound seems not to be competitive with the prospective Katrin bound, the 
neutrino mass bound from SN has a lot merit and importance, and, of course, it should be a source of strong 
investigation. With improvements in the detector efficiency, the number of events per bin would increase and a more 
















6 


competitive bound could be obtained. We remember that efficiency includes several effects, such as, loss of events due 
to light attenuation, fluctuations of the number of photoelectrons, the geometry of the detector etc. These are still 
unknown sources of errors, and their knowledge is crutial for a more detailed work than presented here. We are aware 
that our assumptions are very crude, however, as stated by the LBNE collaboration in m work in understanding the 
physics and choices of detector is still underway. It is necessary to determine, for example, the low energy response 
of the detector, its geometry, its real threshold, the reconstruction of the events, event vertex resolution, the energy 
resolution, absolute event ms time precision and other details. That is why we are being very conservative, but at the 
same time, very optimistic in this work. However, we should stress that considering this time binned analysis, we have 
demonstrated that statistical fluctuation may play a major role in the determination of the neutrino mass bound. We 
point out that a more dedicated MC simulation, with real and proper detector information, with an effective study of 
a possible event-by-event likelihood of detection for this flux of v e of neutronization burst and, finally, the importance 
of the time of the first event detected should be studied deeply in the foreseeable future. 
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